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ABSTRACT 

The origin of the shallow decay segment in the Swift/XRT light curves is still a puzzle. We analyze 
the properties of this segment with a sample of 53 long Swift GRBs detected before Feb., 2007. We 
show that the distributions of its characteristics are log-normal or normal, and its isotropic X-ray 
energy (£7j so ,x) is linearly correlated with the prompt gamma-ray energy, but with a steeper photon 
spectrum except for some X-ray flashes. No significant spectral evolution is observed from this phase 
to the follow-up phase, and the follow-up phase is usually consistent with the external shock models, 
implying that this shallow decay phase is also of external shock origin, likely due to a refreshed external 
shock. Within the refreshed shock model, the data are generally consistent with a roughly constant 
injection luminosity up to the end of this phase tb- A positive correlation between -EjscX and tb also 
favors the energy injection scenario. Among the 13 bursts that have well-sampled optical light curves, 
6 have an optical break around tb and the breaks are consistent with being achromatic. However, 
the other 7 bursts either do not show an optical break or have a break at a different epoch than tb- 
This raises a concern to the energy injection scenario, suggesting that the optical and X-ray emissions 
may not be the same component at least for a fraction of bursts. There are 4 significant outliers in 
the sample, GRBs 060413, 060522, 060607A, and 070110. The shallow decay phase in these bursts 
is immediately followed by a very steep decay after tb, which is inconsistent with any external shock 
model. The optical data of these bursts evolve independently from the X-ray data. These X-ray 
plateaus likely have an internal origin and demand continuous operation of a long-term GRB central 
engine. We conclude that the observed shallow decay phase likely has diverse physical origins. 
Subject headings: radiation mechanisms: non-thermal: gamma-rays: bursts: X-rays 



1. INTRODUCTION 

The observations for the gamma-ray burst (GRB) 
phenomenon with the Swift satellite (Gehrels et al. 
2004) have revolutionized our understanding on this 
phenomenon in many aspects (see recent reviews by 
Meszaros 2006; Zhang 2007). In its first two years of 
operation, the on-board X-ray telescope (XRT; Burrows 
et al. 2004) has accumulated a large set of well-sampled 
X-ray light curves from tens of seconds to days (even 
months) since the GRB triggers. 

The generally accepted GRB models are the relativis- 
tic fireball models (Rees & Meszaros 1994; Meszaros & 
Rees 1997a; Sari et al. 1998; see reviews by Meszaros 
2002; Zhang k, Meszaros 2004; Piran 2005). This model 
invokes a fireball powered by a GRB central engine 
that ejects an intermittent, relativistic outflow. Internal 
shocks from stochastic collisions within the ejecta power 
the observed prompt gamma-rays, and deceleration of 
the fireball by the ambient medium excites a long term 
external forward shock that powers the broad band af- 
terglow (Meszaros & Rees 1997a; Sari et al. 1998). Swift 
data suggest possible late internal shocks that are the 
origin of the erratic late X-ray flares seen in XRT light 
curves (Burrows et al. 2005; Zhang et al. 2006; Dai et al. 
2006; Fan & Wei 2005; King et al. 2005; Proga & Zhang 
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2006; Perna et al. 2006). The XRT light curves gener- 
ally begin with a rapidly decaying segment (Tagliafcrri 
et al. 2005; O'Brien et al. 2006b), which is explained 
as the prompt emission tail due to the so-called "curva- 
ture effect" (Kumar & Panaitescu 2000b; Zhang et al. 
2006; Liang et al. 2006; Yamazaki et al. 2006) 4 . The 
broadband afterglows, which usually decay as a power- 
law with an index of a ~ — 1 (normal decay phase), are 
believed to be related to the external shock. If the exter- 
nal shocks are refreshed by continuous energy injection 
into the blastwave, a shallow decay phase prior of the 
normal decay phase could be observed (Rees & Mezaros 
1998; Dai & Lu 1998a,b; Panaitescu et al. 1998; Sari & 
Meszaros 2000 ; Zhang & Meszaros et al. 2001; Wang & 
Dai 2001; Dai 2004; Granot & Kumar 2006; Panaitescu 
2007; Yu & Dai 2007; see Zhang 2007 for a review). 

As the fireball is decelerated by the ambient medium, 
the normal decay phase transits to a jet-like decay phase 
(with a decay index a ~ —2), when the bulk Lorentz fac- 
tor is degraded to V ~ 6 J 1 , where 9j is the opening angle 
of a conical jet (Rhoads 1997; Sari et al. 1999). There- 
fore, four successive emission episodes are invoked in the 
framework of the fireball models, i.e., prompt gamma- 
ray phase with a tail, shallow decay phase, normal de- 
cay phase, and jet-like decay phase. These power law 
segments, together with erratic X-ray flares, composes 
a canonical X-ray afterglow lightcurve revealed by Swift 
(Zhang et al. 2006; Nousek et al. 2006; O'Brien et al. 

4 Other suggestions to interpret this segment include cooling of 
a hot cocoon surrounding the GRB jet (Pe'er et al. 2006) or a 
highly radiative blast wave (Dermer 2007). 
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2006b). The physical origins of these segments have been 
discussed in the literature (Zhang et al. 2006; Nousek et 
al. 2006; Panaitescu et al. 2006a). Empirically, O'Brien 
et al. (2006a,b) and Willingale et al. (2007) show that 
the data can be fitted by the superposition of a prompt 
emission component and an afterglow component. 

In order to explore the physical origin of this canonical 
afterglow lightcurve, we perform a systematic analysis 
of the Swift XRT data. In the first paper of this series 
(Zhang, Liang, & Zhang 2007, Paper I of this series), 
we have studied the steep decay phase for a sample of 
bright tails, and revealed an apparent hard-to-soft spec- 
tral evolution for some bursts (see also Campana et al. 
2006; Mangano et al. 2007; Butler & Koccvski 2007). 
This paper will focus on the shallow decay phase and the 
subsequent phase. This is motivated by some puzzling 
facts related to the shallow decay phase. For example, si- 
multaneous X-ray/optical observations suggest that the 
break between the shallow and the normal decay seg- 
ments in the X-ray band for some GRBs is chromatic 
(Panaitescu et al. 2006b; Fan & Piran 2006). This is 
inconsistent with the simplest energy injection model. 
One fundamental question is whether X-ray and optical 
afterglows have the same physical origin. Another inter- 
esting fact is that the XRT light curve of GRB 070110 
shows a long-lived plateau followed by an abrupt fall- 
off (the decay slope ~ — 9 with zero time at the trigger 
time). This feature is hard to interpret within the exter- 
nal shock models, and it likely indicates a long-lasting 
central engine emission component (Troja et al. 2007). 

Theoretically, several models have been proposed to in- 
terpret the shallow decay phase (e.g. Zhang 2007 for a re- 
view). Besides the energy injection models (Zhang et al. 
2006; Nousek et al. 2006; Panaitescu et al. 2006a), other 
models include the combination of the GRB tail with 
the delayed onset of the afterglow emission (Kobayashi 
& Zhang 2007); off-beam jet model (Toma et al. 2006; 
Eichler & Granot 2006); pre-cursor model (Ioka et al. 
2006); two-component jet (Granot et al. 2006; Jin et 
al. 2007), varying microphysics parameter model (Ioka 
et al. 2006; Panaitescu et al. 2006b; Fan & Piran 
2006; Granot et al. 2006), etc. The chromaticity of 
some X-ray shallow-to-normal breaks drives several ideas 
that go beyond the traditional external forward shock 
model. For example, Shao & Dai (2007) interpret the X- 
ray lightcurve as due to dust scattering of some prompt 
X-rays, so that it has nothing to do with the external 
shock. Uhm & Beloborodov (2007) and Genet, Daigne 
& Mochkovitch (2007) interpret both X-ray and optical 
afterglow as emission from a long-lived reverse shock. 
Ghiscllini et al. (2007) even suggested that the shallow- 
to-normal transition X-ray afterglows may be produced 
by late internal shocks, and the end of this phase is due 
to the jet effect in the prompt ejecta (see also Nava et 
al. 2007). 

The observational puzzles and theoretical chaos call for 
a systematic understanding of the shallow decay phase 
data for a large sample of GRBs. In particular, it is 
desirable to find out how bad the standard external for- 
ward shock model is when confronted with the data, e.g. 
what fraction of bursts actually call for models beyond 
the standard external forward shock model. This is the 
primary goal of this paper. Data reduction and sam- 
ple selection are presented in §2. The characteristics of 



the shallow decay segment and their relations with the 
prompt gamma-ray phase are explored in §3. In §4, we 
test the external origin of the power law segment fol- 
lowing the shallow decay phase and explore whether or 
not the shallow decay segment is also of external origin. 
Assuming an energy injection model for shallow decay 
phase, we further analyze the the energy injection model 
parameters of these bursts in §5. The relation among 
the isotropic X-ray energy (Ei SOy x), peak energy of the 
prompt gamma-ray vf u spectrum (E p ), and th is inves- 
tigated in §6. The results are summarized in §7 with 
some discussion. Throughout the paper the cosmologi- 
cal parameters H = 71 km s _1 Mpc -1 , Qm = 0.3, and 
Oa = 0.7 have been adopted. 

2. DATA REDUCTION AND SAMPLE SELECTION 

The XRT data are taken from the Swift data archive. 
We developed a script to automatically download and 
maintain all the XRT data. The Heasoft packages, in- 
cluding Xspec, Xselect, Ximage, and Swift data analysis 
tools, are used for the data reduction. We have devel- 
oped an IDL code to automatically process the XRT data 
for a given burst in any user-specified time interval. Our 
procedure is described as follows. The details of our code 
have been presented in Paper I. 

Our code first runs the XRT tool xrtpipeline to repro- 
duce the XRT clean event data, and then makes pile-up 
corrections with the same methods as discussed in Ro- 
mano et al. (2006) (for the Window Timing [WT] mode 
data) and Vaughan et al. (2006) (for the Photon Count- 
ing [PC] mode data). Both the source and background 
regions are annuli (for PC) or rectangular annuli (for 
WT). The inner radius of the (rectangular) annuli are 
dynamically determined by adjusting the inner radius of 
the annuli through fitting the source brightness profiles 
with the King's point source function (for PC) or deter- 
mined by the photon flux using the method described 
in Romano et al 2006 (for WT). If the pipe- up effect is 
not significant, the source regions are in the shape of a 
circle with radius R = 20 pixels (for PC) or of a 40x20 
pixels rectangle (for WT) centered at the bursts' posi- 
tions. The background regions have the same size as the 
source region, but has a distance of 20 pixels away from 
the source regions. The exposure correction is also made 
with an exposure map created by XRT tools xrtexpomap. 
By considering these corrections, the code extracts the 
background-subtracted light curve and spectrum for the 
whole XRT data set. The signal-to-noise ratio is nor- 
mally taken as 3 cr, but it is not rigidly fixed at this value 
and may be flexibly adjusted depending on the source 
brightness. 

With our code we process all the XRT data observed 
between Feb., 2005 and Jan., 2007. We inspect all the 
light curves to identify the beginning (ti) of the shallow 
decay segment and the end (£2) of the decay phase fol- 
lowing the shallow phase (which usually is the normal 
decay phase, but in some cases the decay slope could be 
much steeper). Please note that the selection of t\ and £2 
is guided by eye without a rigid criterion. Generally, t\ 
is taken as the end of the steep decay segment or the be- 
ginning of the observation time, unless significant flares 
or high level emission bumps following the GRB tails 
were observed. The ending time ti is taken as the end of 
the observation time except for GRBs 050416A, 050803, 
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060413, 060908, 060522, 061121,and 070110, which have 
an additional break at later times, and ti is chosen as that 
break time. For example, GRBs 060522 and 070110 have 
a distinct "normal-decay" emission component following 
the sharp decay segment, and i 2 is taken the end of the 
sharp decay. The last data points of GRB 050416A, 
050803, 060413, 060908, and 061121 show a flattening 
feature, which significantly deviates from the power law 
decay trend post tb- We thus do not include those data 
points. 

Physically, the temporal break of an external shock 
origin should be smooth (due to the equal-arrival-time 
effect of a relativistic shell of conical geometry). There- 
fore, a smoothed broken power law is used to fit the light 
curve in the time interval [t\, t 2 ], 

-l/w 



F = F n 



(1) 



t \ 1 / t 

Jb) + \ib, 
where lo describes the sharpness of the break. The larger 
the lo, the sharper the break. In order to constrain w 
it is required that the time interval covers a range from 
ti <C tb to ti 3> tb, and that the light curve around tb 
is well-sampled. The parameter tb is not significantly 
affected by lo, but both a\ and a 2 are. We show the 
comparison of the fitting results with lo = 1 and lo = 3 
for the bursts in our sample (see below) in Fig. 1. We 



find that systematically, ft 
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< at 
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a^- 1 > o^-*. We notice that Willingale et al. (2007) fit 
the XRT light curves with a superposition model of both 
the prompt and afterglow emission components. The de- 
rived a 2 from our fitting with lo = 3 is generally consis- 
tent with their results. We therefore fix lo = 3 in this 
analysis, except for GRBs 060413, 060522, 060607A, and 
070710. The XRT light curves of these bursts abruptly 
drop at tb, and we take lo = 10. We then create a time 
filter array that contains two time intervals of [ii,ib] 
and [tb^] for each burst. By specifying the time fil- 
ter array we run our code again to extract the spectra 
in the two time intervals and derive their photon indices, 
Tx.i and Fx, 2, by fitting the spectra with a power law 
model incorporating with absorptions by both the Milky 
Way Galaxy and the host galaxy, wabs Gal x zwabs host x 
power law (when the redshift is unknown, zwabs host is re- 
placed with the model of wabs). The N^f st value in the 
time-resolved spectral analysis is fixed to the value ob- 
tained from fitting the time-integrated spectrum during 
the whole time span of each burst. 

The tb is roughly considered as the duration of the 
shallow decay phase. As suggested by Lazzati & Begel- 
man (2006) and Kobayashi & Zhang (2007), the zero 
time of the external-origin power-law segments should 
be roughly the BAT trigger time. In our calculation, in 
order to account for the onset of the afterglow we take 
a to as 10 seconds after the GRB trigger. The X-ray 
fluence (Sx) of the shallow decay phase is derived by in- 
tegrating the fitting light curve from 10 seconds post the 
GRB trigger to tb without considering the contributions 
of both early X-ray flares and the GRB tail emissions. 
Since the shallow decay phase has a temporal decay in- 
dex shallower than -1, the results are not sensitive to the 
choice of t . We estimate the uncertainty of Sx with 
a boostrap method based on the errors of the fitting 
parameters, assuming that the errors of the fitting pa- 
rameters, CTi og _F , criogtb, cr ai , and of Gaussian 



distributions. We generate 5 x 10 3 parameter sets of (Fq, 
tb, 011, ol 2 ) from the distributions of these parameters for 
each burst, and then calculate Sx for each parameter set. 
We make a Gaussian fit to the distribution of log Sx and 
derive the central value of logSx and its error <Ji og Sx- 
In our fittings, cui and/or % are fixed for GRBs 050801 
and 060607A. We do not calculate the errors for the two 
bursts. 

We use the following criteria to select our sample. 
First, the XRT light curves have a shallow decay seg- 
ment following the GRB tails. Since the decay slope of 
the "normal" decay phase predicted by the external GRB 
models is generally steeper than 0.75, we require that the 
so-called shallow decay segment has a slope ax.i < 0.75 
at la error. Second, both the shallow decay segment 
and the follow-up segment are bright enough to perform 
spectral analysis. Systematically going through all the 
Swift XRT data before Feb. 2007 we use the above cri- 
teria to compile a sample of 53 bursts. Please note that 
the apparently long GRB 060614 is also included in our 
sample, although it may belong to the short-type bursts 
(Gehrels et al. 2006; Zhang et al. 2007b; Zhang 2006). 
The XRT light curves and the fitting results are shown 
in Fig. 2, and the data are summarized in Table 1. We 
collect the BAT observations of these bursts from GCN 
circular reports, and report them in Table 2. We search 
the optical afterglow data of these bursts from published 
papers and GCN circular reports 5 . We identify a burst 
as optically bright, if three or more detections in the 
UV-optical bands were made. We find that 30 out of the 
53 bursts are optically bright, but only 15 bursts have 
an optical light curve with good temporal coverage. We 
make the Galactic extinction correction and convert the 
observed magnitudes to energy fluxes. We fit these light 
curves with a simple power law or the smooth broken 
power law (lo is also fixed as 3). The fitting results are 
summarized in Table 3. We directly compare the optical 
data with the XRT data in Fig. 2 in order to perform a 
quick visual check of achromaticity of these light curves. 
If multi-wavelength optical light curves are available, we 
show only the one that was observed around the X-ray 
shallow decay phase with the best sampling. Notice that 
the contribution from the host galaxy to the optical light 
curve of GRB 060614 has been removed. 

Twenty-seven out of the 53 GRBs in our sample have 
redshift measurements. Table 4 reports the properties of 
these bursts in the burst rest frame, including the dura- 
tions (T 90 and t h ) and the equivalent-isotropic radiation 
energies (i?i so , 7 and Ei so ,x) m the prompt phase and in 
the shallow decay phase, and the peak energy of the vf v 
spectrum (E p ). The i?i so , 7 and i?i so ,x are calculated by 

InDlS^x) 



iso,(7,X) 



l + z 



(2) 



where S^ is the gamma-ray fluence in the BAT band and 
Sx is the X-ray fluence in the shallow decay phase in 
the XRT band, and Dl is the luminosity distance of the 
source. Due to the narrowness of the BAT band, the 
BAT data cannot well constrain the spectral parameters 
of GRBs (Zhang et al. 2007a). Generally the BAT spec- 
trum can be fitted by a simple power law, and the power 

5 The GRBlog web page (http://grad40.as.utexas.edu/grblog.php) 
has been used. 
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law index T is correlated with E p (Zhang et al. 2007b; 
see also Sakamoto et al. 2007; Cabrera et al.2007) 6 , i.e., 

log E p = (2.76 ± 0.07) - (3.61 ± 0.26) log V. (3) 

We estimate E p with this relation if it is not constrained 
by the BAT data. Wc then calculate the bolometric en- 
ergy Ef so in the 1 — 10 4 keV band with the fc-correction 
method used by Bloom et al. (2001), assuming that the 
photon indices are -1 and -2.3 before and after E p , re- 
spectively (Preece et al. 2000). Both E p and E% are 
listed in Table 4. 

3. THE CHARACTERISTICS OF THE SHALLOW DECAY 

PHASE AND ITS RELATIONS TO THE PROMPT 
GAMMA-RAY PHASE 

We display the distributions of the characteristics 
of the shallow decay phase in Fig. 3. It is 
found that these distributions are consistent with be- 
ing normal/lognormal, i.e. logib/s = 4.09 ± 0.61, 
log S x /erg cm" 2 = -6.52 ± 0.69, r x ,i = 2.09 ± 0.21, 
and ai — 0.35 ±0.35. Quoted errors are at la confidence 
level. 

We investigate the relation of the shallow decay phase 
to the prompt gamma-ray phase. Figure 3 shows t^, 
Sx, Tx,i, and i?i so ,x as functions of T 90 , S 7 , T 7 , and 
Ei SOi7 , respectively. No correlation between T 7 and Tx,i 
is observed. However, Tx.i is larger then T 7 , except for 
some X-ray flashes (XRFs), indicating that the photon 
spectrum of the shallow decay phase is generally steeper 
than that of the prompt gamma-ray phase for typical 
GRBs. It is interesting to note that in contrast to T 7 
Tx,i is narrowly clustered around 2.1 (see also O'Brien 
et al. 2006a), hinting a possible common microscopic 
mechanism during the shallow decay phase. 

From Fig. 3 we find tentative correlations of dura- 
tions, energy fluences, and isotropic energies between 
the gamma-ray and X-ray phases. The best fits yield 
logt b = (0.61 ±0.16) logT 90 + (3.00 ±0.27) (r = 0.48 and 
p = 0.003 for N = 53), logS x = (0.76 ± 0.11) logS^ + 
(-2.33 ± 0.60) (r = 0.70 and p < 10~ 4 for N = 53), 
and \ogE iso . x = (1.00 ± 0.16) log £ iso . 7 + (-0.50 ±8.10) 
(r = 0.79 and p < 10~ 4 for N = 27). It is found 
that tb weakly depends on Tg . However, X-ray fluence 
and isotropic energy are almost linearly correlated with 
gamma-ray fluence and gamma-ray energy, respectively. 
Ei S0 ^ is greater than Ex,iso for most of the bursts, but 
for a few cases Ei SQ ^x is even larger than E^j so . In order 
to reveal possible linear correlations for the quantities in 
the two phases, we define a 2a linear correlation regions 
with y = x + (A ± 2 x a a), where y and x are the two 
quantities in question, and A and a a are the mean and 
its la standard error of the y — x correlation, respec- 
tively. The la regions of the correlations are shown with 
dashed lines in Fig. 4. These results indicate that radi- 
ation during the shallow decay phase is correlated with 
that in the prompt gamma-ray phase. 

4. TESTING THE PHYSICAL ORIGIN OF THE SHALLOW 
DECAY SEGMENT USING THE PROPERTIES OF THE 

FOLLOW-UP SEGMENT 



6 We should point out that this empirical relation is for BAT ob- 
servations only. The origin of this relation is due to the narrowness 
of the BAT instrument. It can be robustly used for those bursts 
whose E p are roughly within the BAT band. 



The leading scenario of the shallow decay phase is a re- 
freshed forward shock due to either a long-term central 
engine or a spread of the ejecta Lorentz factor (Zhang et 
al. 2006; Nousek et al. 2006; Rees & Meszaros 1998; Dai 
& Lu 1998a,b; Zhang & Meszaros 2001; Granot & Kumar 
2006; Yu & Dai 2007). Within such a scenario, the shal- 
low decay phase ends at and transits to a "normal" 
decay phase consistent with the standard external for- 
ward shock models. Three criteria are required to claim 
an energy injection break First, there should be no 
spectral evolution across 4 since energy injection is a 
pure hydrodynamical effect. Second, due to the same 
reason, the break at tf, should be achromatic. Third, the 
power-law decay phase after tf, should comply with the 
standard external shock models. In this section, we test 
whether all three criteria are satisfied with the data. 

Figure 4(a) shows Tx,2 as a function of r x ,i- The 
solid line is Tx,2 = rx,i, and the dashed lines marks 
the 3a region of the equality, which is defined with 
Tx,2 = Txa + [Q ± 3Sg), where Q and Sg are the mean 
and statistical uncertainty (la level) of the difference 
Tx.2 — r Xj i. Please note that Sg does not include the 
observational uncertainty. It statistically describes the 
scatter of Q for the bursts in our sample. We find that 
only one burst GRB 061202 is out of the region. The 
comparison between the distributions of r x .2 and r x ,i is 
shown in Fig. 4(b). Excluding GRB 061202, the two dis- 
tributions are consistent. The Kolmogorov-Smirnov test 
suggests that the significance level of this consistency is 
0.96. These results indicate that r x ,i and Tx,2 for the 
bursts in our sample are globally consistent with each 
other. In order to verify this consistency within observed 
uncertainty for individual bursts, Fig. 4(c) shows the dis- 
tribution of the ratio /i = Q/a, where a 2 — 8T 2 X l +ST X 2 
is the observed uncertainty of Q. A positive value of fj, 
would indicate a hard-to-soft spectral evolution. This 
ratio indicates the significance level of the difference be- 
tween rx,i and Tx,2 for individual bursts within the ob- 
servational uncertainties of the two quantities. As shown 
in Fig. 4(c), most of the bursts (<~ 90%) have \i < 1, 
and only one burst (GRB 061202) has p > 3. These 
results prove that no significant spectral evolution be- 
tween the two phases with a confidence level above 3a is 
observed for the bursts in our sample within the obser- 
vational error, except for GRB 061202. This is consis- 
tent with the expectation of the refreshed shock afterglow 
models. Please note that GRB 061202 shows significant 
hard-to-soft spectral evolution from the shallow to the 
normal decay phases, i.e., from Tx.i = 2.25 ± 0.07 to 
Tx,2 = 3.55 ± 0.44. One caveat for this spectral evo- 
lution is that there is a long observational gap between 
the first epoch in the shallow decay phase (4 x 10 3 to 
2 x 10 4 seconds) and the second epoch in normal decay 
phase (1 x 10 5 ~ 5 x 10 5 seconds) when the spectral in- 
dices are measured. Without detecting the break itself, 
it may be dangerous to draw the conclusion that spec- 
tral variation is clearly seen across tb- We cannot rule 
out the possibility that the plateau extends further and 
drops dramatically before landing onto a normal decay 
segment as is seen in GRBs 060522 and 070110 (see dis- 
cussion below). 

Although the mechanism of energy injection into the 
forward shock could vary (e.g. Rees & Mezaros 1998; Dai 
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& Lu 1998a,b; Zhang & Meszaros et al. 2001; Yu & Dai 
2007), the kinetic energy of the fireball after the energy 
injection is over should be constant and this "normal" de- 
cay phase should be explained with the standard external 
shock models. Without broadband afterglow modeling, 
the "closure relations" between the observed spectral in- 
dex and temporal decay index a present a simple test 
to the models. In Fig. 5, we present ax. 2 as a function 
of spectral index /?x,2 5 where /3x,2 = Tx,2 — 1- The clo- 
sure correlations of the external shock afterglow models 
for different spectral regimes, different cooling schemes, 
different ambient medium properties, and different elec- 
tron distributions (the spectral index p > 2 and p < 2) 
are shown in Fig. 5 (see Table 1 of Zhang & Meszaros 
2004 and reference therein, in particular Sari et al. 1998; 
Chevalier & Li 2000; Dai & Cheng 2001). The fact that 
the observed 0x,2 is greater than 0.5 for the bursts in 
our sample suggests that these X-rays are in the spectral 
regime vx > max(^ m ,i/ c ) (Regime I) or u m < vx < v c 
(Regime II), where v m and v c are the characteristic fre- 
quency and cooling frequency of synchrotron radiation. 
The relation between a and for the spectral Regime I is 
a = (3/3— 1) /2 regardless of the type of the medium (ISM 
or wind medium) . If the X-ray band is in the Regime II, 
we have a = 3/3/2 (for ISM) and a = (3/3 + l)/2 (for 
wind). We define 

D = \a ohs -a(/3 obs )|, 

5 D = ^(6a obs ) 2 + [Sa(0 obs )} 2 , (4) 

where a ohs (Sa° hs ) and a(/3 obs ) are the temporal decay 
slopes (errors) from the observations and from the clo- 
sure relations. The ratio <fi — D/Sd reflects the "near- 
ness" of the data point to the model predictions within 
the error scope. If (f> < 1, we consider that the data 
point goes cross the corresponding closure relation line. 
If 4> < 3, we regard that the model cannot be excluded 
within a 3<r significance level. Those bursts that have 
large uncertainties on both a and may be interpreted 
with more than one models. In this case, we compare 
<j> values derived from these models and take the model 
that gives the smallest (f>. 

As shown in Fig. 5, 24 out of the 53 bursts distribute 
around the line for the spectral Regime I. A group of 
bursts have a decay slope shallower than the model pre- 
diction, but they are slightly below and almost keep 
abreast with the Regime I model line (see also Fig. 5 
of Willingalc et al. 2007). At the 3<r confidence level, 
this model cannot be excluded for these bursts. Eigh- 
teen bursts are consistent with the ISM external shock 
afterglow model in the spectral regime II. 

Six bursts (GRBs 050315, 050318, 050803, 060614, 
051008, and 060906) agree with both the regime I ISM 
jet model and the regime II wind model. The observed 
of these six bursts are ~ 1. The two models are al- 
most degenerate at ~ 1. We therefore use the spectral 
and temporal behaviors of the prior segment to distin- 
guish the two models. Since the observed > 0.5 in 
our sample, the decay slope of the light curves before a 
jet break should be steeper than (3/3)/2 ~ 0.75. From 
Table 1 we can see that the <x\ values are 0.66 ± 0.03, 
0.90 ± 0.23, 0.25 ± 0.03, 0.18 ± 0.06, 0.78 ± 0.11, and 
0.35 ± 0.10, respectively, for the six bursts. So there is 
no confident evidence to claim a jet break within the un- 



certainty of the decay slope for these bursts 7 Since the 
energy injection model the shallow decay slope depends 
on a free parameter q, we tentatively suggest that these 
six cases can be explained with a wind afterglow model 
in the spectral regime II. GRB 060108 is also consistent 
with this model according to our criterion. 

As shown above, the spectral index and temporal decay 
slope of the normal decay phase for most bursts in our 
sample (49 out of 53 bursts) are roughly consistent with 
the closure relations of the external shock models. This 
further favors the idea that the shallow decay segment 
is also of external shock origin, and probably is related 
to a long-term energy injection effect. In this scenario, 
the energy injection break should be achromatic, if the 
multi-wavelength radiations are all from the same emis- 
sion region, presumably the forward shock. We therefore 
inspect the optical light curves of these bursts to examine 
whether the breaks observed in the XRT light curves are 
achromatic. Among these 13 bursts have well-sampled 
optical light curves, as shown in Fig. 2 and Table 3. The 
optical light curves of GRBs 050801, 051109A, 060614, 
060714, 060729, and 061121 show a break around t b , indi- 
cating that the breaks in both the X-ray and the optical 
bands are consistent with being achromatic. However, 
the optical light curves of GRBs 050318, 050319, 050802, 
060124, and 050401 do not have a break around % (see 
also Panaitescu et al. 2006b). They can be well fitted 
by a simple power law model. GRBs 060210 and 060526 
have an optical break, but the breaks are not around tf,. 

GRBs 060413, 060522, 060607A, and 070110 have a 
plateau with a step-like sharp drop (w = 10 is required 
in our data fitting). Except for GRB 060522, the other 
three bursts deviate significantly from any external shock 
afterglow models at 3er significance level. Although the 
sharp drop segment of GRB 060522 is consistent with the 
regime II, wind-jet model, the plateau convincingly rules 
out this model since it cannot be explained as the pre- 
jet segment within the same model. These results sug- 
gest that the sharp drop segment and its prior plateau in 
these bursts are very likely not of external shock origin. 
A direct support to this speculation is that the optical 
lightcurves of these bursts, if available, are all evolve in- 
dependently with respect to the X-ray lightcurves. For 
example, The optical light curve of GRB 060607A rapidly 
increases (with F oc t 3 ) up to a maximum at t = 160 sec- 
onds post the GRB trigger, and then continuously decays 
with an index of -1.18 ±0.02 (Molinari et al. 2007). The 
X-ray lightcurve, on the other hand, shows significant 
flares before 600 seconds, and a plateau lasting from 600 
seconds to ~ 1.2 x 10 4 seconds after the GRB trigger. 
At the end of the plateau, the XRT light curve drops 

7 The possibility that a jet break is temporarily coincident with 
an energy injection break is however not ruled out. The normal 
decay phase could be missed in data fitting, if the normal decay 
segment is short, the data is sparse, or the jet break is not sig- 
nificance. g. Wei & Lu 2000; Kumar & Panaitescu 2000a; Gou 
et al. 2001). One example of this scenario is GRB 060614. Our 
best fit with a smooth broken power law does not reveal a jet-like 
break from a normal decay phase. However, Mangano et al. (2007) 
suggest a normal decay phase between 3.66 X 10 4 and 1.04 X 10 s 
seconds by fitting the light curve with a joint-power-law model 
(the breaks arc guided by eye). They showed that the decay slope 
during the period is 1.03 ± 0.02. This normal decay phase thus 
satisfies a closure relation of the standard forward shock models. 
For a detailed study of jet breaks please see our paper III in the 
series, E.-W. Liang et al. 2007, in preparation. 
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sharply with a 2 — 3.35 ± 0.09. During the plateau in 
the XRT light curve, the optical light curve "normally" 
decays until a significant flare around 2000 seconds. The 
optical light curve is consistent with a external forward 
shock, and the peak is consistent with onset of the after- 
glow (Molinari et al. 2007). The plateau and the sharp 
drop in the XRT light curve of GRB 070110 is similar 
to that of GRB 060607A, but an additional "normal" - 
decay component post the steep fall-off was also observed 
(Troja et al. 2007). The decay slope of this late X-ray 
emission component is similar to that of the optical light 
curve and is likely of external shock origin (sec also GRB 
060522). This reinforces the suggestion that the early X- 
ray plateau is of internal origin and is connected to a 
long-lasting central engine (Troja et al. 2007). A com- 
mon signature of these internal-origin plateaus is that the 
flux almost keeps constant on the plateau but with sig- 
nificant flickering. Although it may not be unreasonable 
to interpret it as late internal shocks (which usually give 
rise to erratic collisions within the ejecta and may power 
X-ray flares), another possibility is that the plateau is 
powered by tapping the spindown energy of the central 
engine, as suggested by Troja et al. (2007). 

5. ENERGY INJECTION BEHAVIOR 

As shown above, the normal decay phase for most of 
the bursts in our sample (49 out of 53) are consistent 
with the external shock models. This suggests that, in 
general, the observed shallow decay phase is also of the 
external origin and may be related to continuous energy 
injection into the fireball. In this section, we assume the 
standard energy injection model and infer from the data 
the parameters of the long-lived central engine. 

We describe the energy injection behavior as L(t) oc 
t~ q (e.g. Zhang & Meszaros 2001) 8 . The difference be- 
tween the decay slopes before and after % depends on the 
observed spectral regime and the type of abient medium, 
which can be summarized as (derived from Table 2 of 
Zhang et al. 2006), 



(p+3)(l-g) 



spectal regime I, ISM and Wind, 
An = I (p+3Ki-g) i S pectal regime II, ISM, 

spectral regime II, Wind, 

where p is the power law index of the electron distribu- 
tion. The p value is derived from the observed spectral 
index, depending on the observed spectral regime. We 
identify the spectral regime for these bursts by comparing 
the observed ax, 2 and ftx,2 with the closure correlations, 
and then derive their q values from Eqs. 5. The distri- 
butions of theses GRBs in the two dimensional q — Aa 
and q — p planes are shown in Fig. 6, along with the 
contours of constant p and Aa lines derived from the 
models (Eq.[5]). No correlation between ot\ and a 2 is 
found. The steepening index Aa is found to vary among 
bursts, with an average of 1.11 ± 0.39. The p values 
range from 2 ~ 3.5 without evidence of clustering (see 

8 Another injection scenario invoking a distribution of the 
Lorcntz factor of the ejecta (Rees & Meszaros 1998) can be ef- 
fectively represented by a long-term central engine (Zhang et al. 
2006). The internal-origin plateaus discussed above suggest that 
at least for some GRBs, a long-lived central engine is indeed in 
operation. 



also Shen et al. 2006). The q values for most bursts are 
around —0.75 to 0.55, with an average of ~ —0.07 ±0.35. 
It is worth commenting that a specific energy injection 
model invoking a spindown pulsar predicts a q value of 
(Dai & Lu 1998a; Zhang & Meszaros 2001). The average 
q value is close to this model prediction. 

6. EMPIRICAL RELATION AMONG E lso - E P - T B 
An empirical relation among Ei so — E p — t' bopt was dis- 
covered with pre-Swift GRBs, where t' b t is the tempo- 
ral break of the optical afterglow light curve in the rest 
frame of the burst (Liang-Zhang relation; Liang & Zhang 
2005) 9 . Willingale et al. (2007) found a similar correla- 
tion to the Ghirlanda relation by assuming that t b in the 
X-ray band are jet breaks. Such a relation may be inter- 
preted as an effective £"i so — E p — t b relation similar to 
the Liang-Zhang relation. In this section we investigate 
the relations among the energies E\ so , i?i so ,x and the 

other two parameters E p and t h . With the data reported 
in Table 4 (for bursts with redshift measurements), wc 
use a multi-variable regression analysis method to search 
for possible dependences of i?i so ,x and Ef so on both E p 

and t h . Our sample is limited to those bursts whose t b 
can be explained as an energy injection break (without 
considering the achromaticity of the break). Among the 
49 bursts 27 have redshift measurements. Since only two 
bursts in the internal-origin plateau sample have redshift 
measurements, we cannot make an analysis to them. Our 
regression model reads 



log Ei so = «o + ki log E p + K 2 log t b 



(6) 



where E p = E p (l + z) and t b = %/(l + z). We measure 
the significance level of the dependences of each variable 
on the model by the probability of a t-test (pt). The 
significance of the global regression is measured by a F- 
test (with a chance probability pf)- Statistically, a ro- 
bust statistical analysis requires the chance probability 
to be less than 10~ 4 . Our multiple regression analysis to 
E iso ^(E p ,t h ) shows «o = 44.0 ± 1.1 (with p t < 10- 4 ), 
m = 1.82 ± 0.33 (with p t < 10~ 4 ), and k 2 = 0.61 ± 0.18 
(with p t = 3 x 10~ 3 ). The p F is < 10~ 4 . These results 
suggest a strong correlation between -Ei so ,x and E p and a 
tentative correlation between i?i so ,x and t h . On the other 
hand, our multiple regression analysis to E^ so 1 {E p , t b ) 
gives k = 48.3 ± 0.8 (with p t < 10~ 4 ), «i = 1.70 ± 0.25 
(withp t < 10~ 4 ), and k 2 = 0.07±0.13 (with p t = 0.486). 
The pf is < 10~ 4 . It is found that the correlation be- 
tween Ef sol and E p is significant, but statistically no 
correlation between Ef so and t h can be claimed. 

With the relation discovered by Willingale et al.(2007), 
one would naively expect a multi-variable correlation 
among Ef so — E p — t b . According to our results, a 
significant correlation among Ef so — E p — t b could be 
indeed claimed with a chance probability pp < 10~ 4 . 

9 If if, opt is interpreted as a jet break, then the relation is sim- 
ilar to the Ghirlanda-rclation (Ghirlanda et al. 2004). However, 
the more empirical Liang-Zhang relation allows more freedom to 
understand the origin of the breaks. 
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However, this correlation is dominated by the correla- 
tion of E^ so - E p only (with a p t < 10~ 4 ), which is 
essentially the Amati- relation (Amati et al. 2002). The 
p t of the dependence between E% and t h is 0.486. This 
strongly rules out such a dependence. Therefore, we 
suspect that the apparent relation found by Willingale 
et al. (2007) would be intrinsically a manifestation of 
the Amati-relation. A similar conclusion has been also 
achieved by Nava et al. (2007). The tb essentially did 
not enter the problem, since the distribution of tb is nar- 
rower than that of tb opt as discussed in Liang & Zhang 
(2005). 

It is interesting to note the dependence -Ei so ,x oc 
£ i o.6i±o.i8 ^his j g m g nar p contrast to the Liang-Zhang 
relation, in which -Ei so , 7 oc t h ~^ 4 was discovered. In 

order to compare the i?i so ,x — E — t h correlation with 
the Liang-Zhang relation in a 2-dimcnsional plane, we 
define S = log_E; so — K 2 logi b , and show Ex and £ 7 
as a function of log.Ep in Fig. 8. We observe that the 

Eiso.x — E p — t h correlation is significant, but it has a 
larger scatter than the Liang-Zhang relation. Although 
we can not rule out the possibility that the large disper- 
sion is intrinsic, the observational uncertainties of both 
Ex,iso and E p could make such a dispersion. Figure 8 
evidently shows that the Ei SOt x—E p —tb correlation is dif- 
ferent from the Liang-Zhang relation. This suggests that 
% and tb.opt may have distinct physical origins. The pos- 
itive correlation between i?i so ,x and % is consistent with 
energy injection origin of %, namely, a longer injection 
episode gives more energy. The negative correlation be- 
tween £7i so , 7 and ib.opt may step back to the standard 
energy reservoir argument of Frail et al. (2001), which 
suggests a connection between ib.opt and the jet opening 
angle of the outflow. 

7. CONCLUSIONS AND DISCUSSION 

We have presented a comprehensive analysis of the 
Swift XRT light curves of long GRBs, focusing on the 
properties of the shallow decay phase and its relation 
with the follow-up decay phase. Our sample includes 53 
bursts whose X-ray emissions are bright enough to per- 
form spectral and temporal analyses for both phases. We 
summarize our results as follows. 

(1) We find that the distributions of the characteristic 
properties of the shallow decay phase are log-normal or 
normal, i.e., logib/s = 4.09 ± 0.61, log .Sx /erg cm~ 2 = 
-6.52 ± 0.69, r x ,i = 2.09 ± 0.21, and ai = 0.35 ± 0.35 
(quoted errors are at la confidence level). 

(2) The E iso _x of the shallow decay phase is lin- 
early correlated with the prompt gamma-ray phase, i.e., 
\ogE x , tso = (1.00 ± 0.16) log E iso>1 - (0.5 ± 8.12) (with 
a Spearman correlation coefficient r = 0.79 and a chance 
probability p < 10~ 4 ). The spectrum of the shallow de- 
cay phase is softer than the prompt gamma-ray phases, 
except for some typical XRFs. 

(3) Except for GRB 061202, no spectral evolution is 
observed during the transition from the shallow decay 
to the follow-up decay phases. The post break phase in 
most bursts is consistent with the closure relations of the 
external shock models. Six out of the 13 bursts with well- 
sampled optical light curves show an achromatic break in 
both X-ray and optical bands, but the other 7 cases either 



do not show any break or have a break at a different 
epoch in the optical band. This poses an issue to explain 
tb of these bursts as the end of the energy injection phase. 

(4) Four bursts (GRBs 060413, 060522, 060607A, and 
070110) in our sample have a post-break phase signifi- 
cantly deviate from the external shock models. The de- 
cay indices are much steeper than model requirements. 
The optical light curves of the latter two bursts evolve 
distinctly from the X-ray light curves. We suggest that 
the X-ray and optical emissions of these bursts are from 
different emission sites, and the X-ray plateaus are of in- 
ternal origin and demand a long-live emission component 
from the central engine. 

(5) Within the scenario of the refreshed external 
shocks, the average energy injection index q ~ 0, sug- 
gesting a roughly constant injection luminosity from the 
central engine. 

(6) With a sub-sample of 27 bursts with known red- 
shifts that satisfy the closure relations of the standard 
external fireball models, we discover an empirical multi- 
variable relation among -Ei SOj x, E , and t' h (Eq.[7]), which 
is distinctly different from the i?i so , 7 — E p — t' b t relation 
discussed by Liang & Zhang (2005). 

(7) There is no significant correlation between t' b and 

the other parameters E- lson and E p (unlike t' hopt ). This 
suggests that the apparent £j j7 — E p relation by assum- 
ing a jet origin of tb (Willingale et al. 2007) is likely a 
manifestation of the Amati-relation. 

These results suggest that the shallow decay segment 
observed in most bursts is consistent with having an ex- 
ternal forward shock origin, probably due to a continuous 
energy injection into the forward shock from a long-lived 
central engine. Therefore, the scenarios that completely 
abandon the external shock models (e.g. Ghiscllini ct 
al. 2007; Genet et al. 2007; Uhm & Beloborodov 2007; 
Shao & Dai 2007) may not be demanded by the data, 
and these models need to explain the apparent consis- 
tency of the a — (3 data with the simple closure relations 
of the forward shock models. 

Since the energy injection break is due to a hydrody- 
namic effect, achromatism is one of a key feature of the 
model. Although about half cases satisfy such a con- 
straint, at least some X-ray breaks are chromatic (even 
if the post break segment is well consistent with the stan- 
dard afterglow model) . This poses a great issue to argue 
that these X-ray breaks are energy injection breaks. In- 
voking different emission regions (e.g. Zhang & Meszaros 
2002) may solve the problem, although more detailed 
modelling is needed. Crossing of a cooling break would 
also result in a temporal break, but it would also lead 
to a change of the spectral index by ~ 0.5. From Table 
1 we find that the changes of the X-ray spectral indices 
across the break of these bursts are 0.01 ±0.10 (050318), 
0.04±0.10 (050319), 0.08±0. 12(050401), and 0.03±0.09 
(050802). These results confidently rule out such a pos- 
sibility. Genet et al. (2007) account for these chromatic 
breaks as due to a long-live reverse shock in which only a 
small fraction of electrons are accelerated. The big issue 
of such an interpretation is how to "hide" the emission 
from the forward shock, which carries most of the energy. 

Assuming a simple energy injection law L(t) oc t~ q , 
we find that averagely speaking the injection luminosity 
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could be almost a constant. This places some constraints 
on the physical models of the energy injection models. 
The constant injection luminosity agree with the expec- 
tation of the energy injection model from a central pulsar 
(Dai & Lu 1998a; Zhang & Meszaros 2001), suggesting 
that the pulsar injection model may be consistent with 
the data at least for some GRBs (see also Grupe et al. 
2007; Fan & Xu 2006; De Pasquale et al. 2007; Yu et al. 
2007). 

The temporal decay slopes of some bursts following the 
shallow decay phase are shallower than the model pre- 
dictions [Fig. 5, see also Fig. 5(a) in Willingale et al. 
2007]. This discrepancy may be alleviated by different 
ways. First, as shown in Fig. 1, a-i could be systemati- 
cally steeper if a smoother broken power law model (with 
smaller lu) is adopted. Second, theoretically the temporal 
breaks involving external shocks are usually not sharp. 
Other effects, such as delay of transfer of the fireball en- 
ergy to the forward shock (Kobayashi & Zhang 2007), 
and the structured jet effect (Zhang et al. 2004; Kumar 
et al. 2004; Yamazaki et al. 2006) would modify the 
simplest closure relations to make a band rather than a 
line in the a — (3 plane. 

One interesting conclusion from this study is that at 
least for a small fraction of bursts (e.g. GRBs 060413, 
060522, 060607A, and 070110), the observed shallow de- 
cay phase is likely of internal origin. This is another 
component other than X-ray flares that are possibly of 
internal origin. Contrary to the erratic X-ray flares, this 
component has a smoother light curve with flickering, 
likely due to a steady component from the central en- 
gine. A possible energy source for such a component 
would be the spin energy from the central engine, and 
an internal dissipation of the spindown power may be 
the origin (e.g. Troja et al. 2007). Tapping of the 



rotation energy is likely through magnetic fields, either 
through dipolar spindown for a central millisecond pul- 
sar (Usov 1992; Dai & Lu 1998a,b; Zhang & Meszaros 
2001) or by the Blandford-Znajek mechanism for a black 
hole central engine (Blandford & Znajek 1977; Meszaros 
& Rees 1997b; Li 2000). If one accepts that such a com- 
ponent is common among bursts, one can speculate that 
the observed early X-ray emission is the sum of differ- 
ent emission components. The competition among these 
components shapes the variety of the X-ray light curves 
one observes. Depending on the relative importance of 
the internal and external components, the shallow de- 
cay segment could be possibly dominated by either the 
radiation from the refreshed shocks or by the steady ra- 
diation component from the internal dissipation of the 
central engine. In the former scenario, the shallow decay 
phase transits to a normal decay phase that is consistent 
with the external shock models. In the later scenario, 
the emission level of the underlying afterglow compo- 
nent is weaker than that from the emission component 
of the central engine, so that one needs a steep dropoff 
from the plateau to land onto the external shock emission 
component. 
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The starting and ending time of our lightcurve fitting 

b 9 

The break time and the decay slopes before and after the break, and the fitting x (degrees of freedom). 

C Thc X-ray flucncc (in units of 10 — 7 erg cm -2 ) of the shallow decay phase calculated by integrating the fitting light curve from 10 seconds post 

the GRB trigger to t h . 

d 



The X-ray photon indices before and after t^. 
The fitting results of these bursts have an unaccepted reduced % 2 due to significant flicking. 
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Shallow Decay Segment of the Swift/XKT Data 



TABLE 2 

BAT OBSERVATIONS AND REDSHIFTS OF OUR SAMPLE 



GRB 


Ts)o(s) 


S-/ 3 - 


7^ b 


r> b 


E p (keV)° 


rcf B AT c 


2 


ref z c 


050128 


13.8(2.0) 


45.00(5.00) 


1.50(0.05) 




133.1(30.2) 


2992 






050315 


96.0(10.0) 


28.00(3.00) 


1.28(0.00) 


2.20 


37.0(8.0) 


3105 


1.95 


3101 


050318 


32.0(2.0) 


21.00(2.00) 


2.10(0.11) 




39.5(12.4) 


3134 


1.44 


3122 


050319 


10.0(2.0) 


8.00(0.80) 


1.25(0.00) 


2.15 


28.0(6.0) 


3119, Z07 


3.24 


3136 


050401 


33.0(2.0) 


140.00(14.00) 


1.15(0.00) 


2.65 


132.0(16.0) 


3173, Z07 


2.90 


3176 


050416A 


2.4(0.2) 


3.80(0.40) 


1.00(0.00) 


3.22 


16.0(3.0) 


3273 


0.65 


3542 


050505 


60.0(2.0) 


41.00(4.00) 


1.50(0.10) 




133.1(41.0) 


3364, Z07 


4.27 


3368 


050713B 


75.0(7.5) 


82.00(10.00) 


1.00(0.13) 




109.0(32.0) 


3600, Z07 






050726 


30.0(3.0) 


43.00(7.00) 


1.00(0.16) 




984.0(200.0) 


3682, Z07 






050801 


20.0(3.0) 


4.40(1.00) 


1.40(0.00) 


2.00 


33.0(7.0) 


3730, Z07 






050802 


13.0(2.0) 


28.00(3.00) 


1.12(0.00) 


2.48 


118.0(77.0) 


3737, Z07 


1.71 


3749 


050803 


110.0(11.0) 


39.00(3.00) 


1.05(0.10) 




150.0(68.0) 


3757, Z07 


0.42 


3758 


050822 


102.0(2.0) 


34.00(3.00) 


1.00(0.00) 


2.48 


36.0(7.0) 


3856 
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280 0f28 01 


540 oono 001 


98f0 091 




865 0( 1 78 0) 


4077 Z07 






051016B 


4.0(0.1) 


1.70(0.20) 


2.38(0.23) 




25.2(11.2) 


4104 


0.94 


4186 


051109A 


36.0(2.0) 


21.00(3.00) 


1.50(0.20) 




133.1(69.0) 


4217 


2.35 


4221 


060105 


55.0(5.0) 


182.00(4.00) 


1.11(0.03) 




394.8(75.2) 


4435 






060108 


14.4(1.0) 


3.70(0.40) 


2.01(0.17) 




46.3(18.1) 


4445 


2.03 


4539 


060109 


116.0(3.0) 


6.40(1.00) 


1.96(0.25) 




50.7(26.3) 


4476 






060124 


800.0(10.0) 


110.00(10.00) 


1.17(0.27) 




326.5(277.3) 


4601 


2.30 


4592 


060204B 


134.0(5.0) 


30.00(2.00) 


0.82(0.40) 




96.8(41.0) 


4671 






060210 


255.0(10.0) 


77.00(4.00) 


1.52(0.09) 




126.9(36.7) 


4748 


3.91 


4729 


060306 


61.0(2.0) 


22.00(1.00) 


1.85(0.10) 




62.4(18.7) 


4851 






060323 


18.0(2.0) 


5.70(0.60) 


1.53(0.17) 




124.0(55.3) 


4912 






060413 


150.0(10.0) 


36.00(1.00) 


1.67(0.08) 




90.4(24.5) 


4961 






060428A 


39.4(2.0) 


14.00(1.00) 


2.04(0.11) 




43.9(13.8) 


5022 






060502A 


33.0(5.0) 


22.00(1.00) 


1.43(0.08) 




158.2(43.4) 


5053 


1.51 


5052 


060507 


185.0(5.0) 


45.00(2.00) 


1.83(0.10) 




64.9(19.3) 


5091 






060510A 


21.0(3.0) 


98.00(5.00) 


1.55(0.10) 




118.3(36.1) 


5108 






060522 


69 0f5 01 


n oon 001 


1 59f0 151 




107 9?42 71 

-L\J t -CM^^. 1 J 


5153 


5.11 


5155 


060526 


13.8(2.0) 


4.90(0.60) 


1.66(0.20) 




92.3(44.5) 


5174 


3.21 


5170 


060604 


10.0(3.0) 


1.30(0.30) 


1.90(0.41) 




56.7(46.1) 


5214 


2.68 


5218 


060607A 


100.0(5.0) 


26.00(1.00) 


1.45(0.07) 




150.5(38.6) 


5242 


3.08 


5237 


060614 

WV7Ui^7 J. i 


102 Of 5 01 


21 7 00(A 00) 


2 1 3f0 04) 




37 5f9 91 


5256 


0.125 


5275 


060707 


68.0(5.0) 


17.00(2.00) 


0.66(0.63) 




66.0(25.0) 


5289 


3.43 


5298 


060708 


9.8(1.0) 


5.00(0.40) 


1.68(0.12) 




88.4(29.4) 


5395 






060714 


115.0(5.0) 


30.00(2.00) 


1.99(0.10) 




48.0(14.5) 


5334 


2.71 


J06 


060719 


55.0(5.0) 


16.00(1.00) 


2.00(0.11) 




47.1(14.7) 


5349 






060729 


116.0(10.0) 


27.00(2.00) 


1.86(0.14) 




61.2(21.7) 


5370 


0.54 


5373 


060804 


16.0(2.0) 


5.10(0.90) 


1.78(0.28) 




71.8(43.7) 


5395 






060805A 


5.4(0.5) 


0.74(0.20) 


2.23(0.42) 




31.8(23.2) 


5421 






060807 


34.0(4.0) 


7.30(0.90) 


1.57(0.20) 




112.9(56.6) 


5403 






060813 


14.9(0.5) 


55.00(1.00) 


0.53(0.15) 




192.0(19.0) 


5443,5446 






060814 


146.0(10.0) 


150.00(2.00) 


1.56(0.03) 




115.6(24.3) 


5459 






060906 


43.6(1.0) 


22.10(1.40) 


2.02(0.11) 




45.5(14.2) 


5538 


3.68 


5535 


060908 


19.3(0.3) 


29.00(1.00) 


1.33(0.07) 




205.5(53.4) 


5551 


2.43 


5555 


060912 


5.0(0.5) 


13.00(1.00) 


1.74(0.09) 




77.9(22.3) 


5561 


0.94 


5617 


061021 


46.0(1.0) 


30.00(1.00) 


1.31(0.06) 




217.1(52.4) 


5744 






061121 


81.0(5.0) 


137.00(2.00) 


1.41(0.03) 




166.5(33.2) 


5831 


1.31 


5826 


061202 


91.0(5.0) 


35.00(1.00) 


1.63(0.07) 




98.6(25.4) 


5887 






061222A 


72.0(3.0) 


83.00(2.00) 


1.39(0.04) 




175.3(36.8) 


5964 






070110 


85.0(5.0) 


16.00(1.00) 


1.57(0.12) 




112.9(38.4) 


6007 


2.35 


6010 


070129 


460.0(20.0) 


31.00(3.00) 


2.05(0.16) 




43.1(16.1) 


6058 







REFERENCES. — J06: Jakobsson et al. 2006(c); Z07: Zhang et al. 2007a; 2992: Cummings et al. (2005) : 3101: Kelson and Berger 

(2005) ; 3105: Krimm ct al. (2005c); 3119: Krimm ct al. (2005a); 3122: Berger and Mulchaey (2005); 3134: Krimm ct al. (2005b); 
3136: Fynbo ct al. (2005a); 3173: Sakamoto ct al. (2005a); 3176: Fynbo ct al. (2005c); 3273: Sakamoto ct al. (2005c): 3364: 
Hullinger et al. (2005a); 3368: Berger ct al. (2005b); 3542: Ccnko ct al. (2005): 3600: Parsons ct al. (2005a); 3682: Barthelmy et 
al. (2005a); 3730: Sakamoto ct al. (2005b); 3737: Palmer ct al. (2005); 3749: Fynbo ct al. (2005b): 3757: Parsons (2005b); 3758: 
Bloom ct al. (2005); 3856: Hullinger ct al. (2006b); 4077: Barthelmy ct al. (2005b); 4104: Barbicr ct al. (2005); 4186: Soderberg 
ct al. (2005); 4217: Fcnimorc et al. (2005): 4221: Quimby ct al. (2005); 4435: Markwardt et al. (2006a); 4445: Sakamoto et al. 
(2006b); 4476: Palmer ct al. (2006a): 4539: Melandri ct al. (2006); 4592: Ccnko ct al. (2006a); 4601: Lamb ct al. (2006); 4671: 
Markwardt ct al. (2006c); 4729: Cucchiara ct al. (2006a); 4748: Sakamoto et al. (2006a); 4851: Hullinger ct al. (2006); 4912: Parsons 
ct al. (2006b); 4961: Barbicr ct al. (2006a): 5022: Markwardt ct al. (2006b); 5052: Cucchiara ct al. (2006b); 5053: Parsons ct al. 
(2006c): 5091:Barbier et al. (2006b): 5108: Barbicr ct al. (2006c); 5153: Krimm ct al. (2006c); 5155: Ccnko ct al. (2006a); 5162: 
Berger and Gladdcrs (2006a); 5170:Berger and Gladders (2006b); 5174: Markwardt et al. (2006d); 5214: Parsons ct al.(2006e); 5218: 
Castro-Tirado ct al. (2006); 5237: Lcdoux ct al. (2006); 5242:Tuellcr ct al. (2006b); 5256: Barthelmy, S. ct al. (2006a); 5275: Price 
et al.(2006): 5289: Stamatikos ct al. (2006b); 5295: Fenimorc ct al. (2006a); 5298: Jakobsson ct al. (2006a); 5320: Jakobsson et 
al. (2006c): 5334: Krimm ct al. (2006d); 5349: Sakamoto ct al. (2006c); 5370: Parsons ct al. (2006a); 5373: Thocnc ct al. (2006); 
5395: Tucllcr et al. (2006c); 5403: Barbicr ct al. (2006d); 5421: Barthelmy ct al. (2006b): 5443: Cummings ct al. (2006); 5446: 
Golenetskii ct al. (2006); 5459: Stamatikos ct al. (2006a); 5535: Vrccswijk (2006); 5538: Sato ct al. (2006); 5551: Palmer ct al. 
(2006c); 5555: Rol et al. (2006); 5561: Parsons et al. (2006d); 5617: Jakobsson (2006d); 5744: Palmer ct al. (2006b); 5826: Bloom 

(2006) ; 5831:Fenimorc ct al. (2006b): 5887: Sakamoto ct al. (2006d); 5964: Tucllcr ct al. (2006a); 6007: Cummings ct al. (2007); 
6010: Jaunscn ct al. (2007); 6058: Krimm ct al. (2006c). 

a The observed gamma-ray fluencc and its error in 15-150 kcV band, in units of 10 — ^ erg cm - ^. 

The spectrum of the prompt gamma-rays is generally fitted by a simple power-law, f u oc v ^ 1 ^ . Wc estimate the E p for these 
bursts with the empirical relation between T-^ and Ep for the BAT observations. Few cases are fitted with a cut-off power-law or 
the Band function (the photon indices before and after the break energy is ^-y^l and r<y,2) respectively). 
C Rcferences for BAT data and redshift data. 
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TABLE 3 

The optical observations and our fitting results 



GRB 


tl (ks) 


<2 (ks) a 


*6,o(ks) 


b 

«o,i 


b 


X /(dot) 


ret 


05031 8 


3.23 


22.83 




84(0 22") 




5(11 


(1) 


050319 


2.00 


204.74 




0.42(0.02) 




11(16") 




050401 


0.06 


1231.18 




0.80(0.01) 




43(12) 


(5)-(7) 


050801 


0.02 


9.49 


0.19(0.02) 


-0.02(0.07) 


1.10(0.02) 


26(42) 


(8) 


050802 


0.34 


127.68 




0.85(0.02) 




50(10) 


(9)-(ll) 


051109A 


0.04 


20170.00 


21.80(10.95) 


0.66(0.02) 


1.10(0.08) 


106(42) 


(12) 


060124 


3.34 


1979.30 




0.85(0.02) 




11(19) 


(13)-(14) 


060210 


0.09 


7.19 


0.70(0.18) 


0.01(0.24) 


1.23(0.08) 


5(12) 


(15)-(16) 


060526 


0.06 


893.55 


84.45(5.88) 


0.67(0.02) 


1.80(0.04) 


116(56) 


(17) 


060607A 


0.07 


13.73 


0.16(fixed) 


-3.07(0.25) 


1.18(0.02) 


92(35) 


(18) 


060614 


1.54 


934.36 


39.09(1.71) 


-0.40(0.05) 


2.16(0.03) 


114(24) 


(19)-(21) 


060714 


3.86 


285.87 


1.00 


O.Ol(fixed) 


1.41(0.03) 


35(11) 


(22)-(26) 


060729 


20.00 


662.39 


43.29(5.15) 


-0.37(0.34) 


1.34(0.06) 


36(27) 


(27) 


061121 


0.26 


334.65 


1.70(0.73) 


0.17(fixed) 


0.99(0.05) 


18(23) 


(28) 


070110 


0.66 


34.76 




0.43(0.08) 




1(4) 


29 



References. — (l)Still et al. (2005); (2) Quimby et al. (2006); (3) Huang et al. (2007); (4) Mason et al. 
(2006); (5) Dc Pasquale (2006); (6) Rykoff et al. (2005); (7) Watson ct al. (2006); (8) Rykoff ct al. (2006); 
(9) McGowan ct al. (2005a); (lO)Pavlonko ct al.(2005); (11: McGowan ct al . (2005b) ;( 12) Yost ct al. (2007); 
(13) Misra ct al. (2007); (14): Romano ct al. 2006; (15) Curran ct al. (2007) ; (16) Stanek et al.(2007); (17) 
Dai ct al. 2007; (IS) Molinari ct al.(2007); (19)Fynbo ct al.(2006); {20)Della Vallc ct al.(2006); (21)Gal-Yam 
ct al.(2006); (22) Asfandyarov ct al.(2006); (23)Rumyantscv ct al.(2006); (24) Pavlenko ct al. (2006a); (25) 
Jakobsson et al. (2006c); {26)Pavlenko et al. (2006b); (27) Grupc ct al.(2007); (28) Page et al.(2007); (29) Troja 
ct al.(2007). 

a Thc time interval concerned in our fitting 

^For a smooth, broken power law fit, q, ctQ \ . ct0 2 arc the break time and the decay slopes before and 
after the break. For a simple power law fit, the decay index and its error arc shown in column olq ]_ . In order 
to make the fittings more reasonable, we assume an error of 0.1 mag for those data points without observational 
error available. 



TABLE 4 

Rest-frame properties of the bursts with known redshifts in our sample 



GRB 


l°g ^(S) 


io g i?;(kcv) 


log_E ia0i7 (erg) 


lo g s f so , 7 ( er g) 


logE is0i x(erg) 


lo s* b ( s ) 


050315 


1.51(0.05) 


2.0(0.1) 


52.41(0.05) 


52.84 


51.94(0.10) 


4.83(0.07) 


050318 


1.12(0.03) 


2.0(0.1) 


52.04(0.04) 


52.38 


51.88(0.46) 


4.03(0.20) 


050319 


0.37(0.09) 


2.1(0.1) 


52.24(0.04) 


52.69 


52.03(0.49) 


4.01(0.51) 


050401 


0.93(0.03) 


2.7(0.1) 


53.41(0.04) 


53.69 


52.82(0.06) 


3.77(0.06) 


050416A 


0.16(0.04) 


1.4(0.1) 


50.62(0.05) 


51.00 


49.69(0.26) 


2.88(0.28) 


050505 


1.06(0.01) 


2.8(0.1) 


53.14(0.04) 


53.51 


52.62(0.13) 


3.90(0.09) 


050802 


0.68(0.07) 


2.5(0.3) 


52.31(0.05) 


52.62 


51.86(0.11) 


3.61(0.06) 


050803 


1.89(0.04) 


2.3(0.2) 


51.24(0.03) 


51.67 


50.29(0.04) 


3.85(0.03) 


051016B 


0.32(0.01) 


1.7(0.2) 


50.59(0.05) 


50.95 


50.99(0.22) 


4.82(0.15) 


051109A 


1.03(0.02) 


2.6(0.2) 


52.43(0.06) 


52.82 


52.65(0.19) 


4.44(0.13) 


060108 


0.68(0.03) 


2.1(0.2) 


51.56(0.05) 


51.89 


51.20(0.14) 


4.34(0.15) 


060124 


2.38(0.01) 


3.0(0.4) 


53.13(0.04) 


53.72 


53.08(0.18) 


4.72(0.09) 


060210 


1.72(0.02) 


2.8(0.1) 


53.36(0.02) 


53.72 


53.18(0.12) 


4.38(0.09) 


060502A 


1.12(0.07) 


2.6(0.1) 


52.10(0.02) 


52.54 


51.82(0.10) 


4.81(0.09) 


060522 


1.05(0.03) 


2.8(0.2) 


52.69(0.04) 


53.03 


51.29(0.19) 


2.16(0.05) 


060526 


0.52(0.06) 


2.6(0.2) 


52.02(0.05) 


52.36 


51.21(0.26) 


3.38(0.25) 


060604 


0.43(0.13) 


2.3(0.4) 


51.32(0.10) 


51.64 


50.88(0.37) 


3.27(0.26) 


060607A 


1.39(0.02) 


2.8(0.1) 


52.72(0.02) 


53.12 


52.84(0.01) 


4.09(0.01) 


060614 


1.96(0.02) 


1.6(0.1) 


50.90(0.01) 


51.24 


49.25(0.05) 


4.70(0.03) 


060707 


1.19(0.03) 


2.5(0.2) 


52.61(0.05) 


52.98 


51.44(1.36) 


3.95(1.06) 


060714 


1.49(0.02) 


2.3(0.1) 


52.69(0.03) 


53.01 


51.95(0.14) 


3.57(0.11) 


060729 


1.88(0.04) 


2.0(0.2) 


51.31(0.03) 


51.65 


51.35(0.02) 


4.86(0.02) 


060906 


0.97(0.01) 


2.3(0.1) 


52.78(0.03) 


53.09 


51.57(0.13) 


3.62(0.10) 


060908 


0.75(0.01) 


2.8(0.1) 


52.59(0.01) 


53.07 


51.12(0.21) 


2.32(0.16) 


060912 


0.41(0.04) 


2.2(0.1) 


51.47(0.03) 


51.83 


50.21(0.17) 


3.05(0.12) 


061121 


1.54(0.03) 


2.6(0.1) 


52.78(0.01) 


53.23 


51.83(0.13) 


3.90(0.08) 


070110 


1.40(0.03) 


2.6(0.1) 


52.31(0.03) 


52.68 


52.19(0.03) 


4.31(0.01) 
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Fig. 1. — Comparison of the fitting results of smooth broken power law models with different sharpness parameters: ui = 1 and ui = 3. 
The solid lines stand for the equality of the two quantities. 
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Fig. 2 (continued) 

Fig. 2. — The XRT light curves (dots) for the bursts in our sample. The solid lines are the best fits with the smooth broken power law 
for the shallow decay phase and its follow-up decay phase (usually the "normal" decay phase). The fitting x 2 an d degrees of freedom are 
shown in each plot. The optical light curves are shown by opened triangles, if they are available. They are fitted by a smooth broken power 
law or a simple power law as displayed by dashed lines. 
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Fig. 3. — Distributions of the characteristics of the shallow decay segment for the bursts in our sample. The dashed lines are the fitting 
results with Gaussian functions. 
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Fig. 4. — The correlations between the data of the shallow decay phase and the prompt gamma-ray phase. The solid line in each panel 
is the best fit. The dashed lines mark a 2a region defined as y = x + (A ± 2 X a a), where y and x are the quantities in the y and rr-axes, 
respectively, and A and a a are the mean and its la standard error of y — x, respectively. The dash-dotted line is y = x. 




Fig. 5. — The comparison between Tx.i and Tx,2: (a) Tx,i vs. Fx, 2- The solid lines marks the equality line Tx,2 = Tx,i and its 2a 
region; (b) The histograms of Tx,i and Tx,2; (c) The distribution of /i. 
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Fig. 6. — The temporal decay index ax, 2 as a function of the spectral index /?x 2 for the post-break segment as compared with the 
closure correlations of various external shock afterglow models: (1) v > max(i/ c , u m )\ (2) u m < u < v c (ISM, slow cooling); (3) v m < v < v c 
(Wind, slow cooling) (4) v > v c (Jet, slow cooling) (5) v m < v < v c (Jet, slow cooling). The solid lines are those for electron distribution 
index p > 2, and the dashed lines are for p < 2. The so/id dots represent the bursts whose ox, 2 and /3x,2 satisfy the models (1) and (2), 
and the open dots represent those bursts can be explained with the model (3). The stars are those bursts that significantly deviate from 
the external shock afterglow models including 060522 (see discussion in the text). 
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Fig. 7. — The distributions of the 49 GRBs whose normal decay phases are consistent with the external shock models in the q — Act 
(panel a) and q — p (panel b) planes along with the model predictions for v > max(i/ c , v m ) (solid lines) and v m < v < v c (dashed lines). 
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Fig. 8. — Comparison of the Ej SOj x — E p — t h relation with the bursts in our sample (the solid dots; the solid and dashed lines mark 
the best fit in 3(7 level) with the Liang-Zhang relation derived with pre-Swift bursts (the open triangles and the dotted line is the best fit; 
Liang & Zhang 2005), where T,x = logi?i SOi x — K 2* b and S 7 = logi?i SOi7 — K2* b q- 



